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Abstract 

A measurement of the underlying event (UE) activity in proton-proton collisions at a 
centre-of-mass energy of 7 TeV is performed using Drell-Yan events in a data sample 
corresponding to an integrated luminosity of 2.2 fb^^, collected by the CMS experi- 
ment at the LHC. The activity measured in the muonic final state (qq — > ji^ji^) is 
corrected to the particle level and compared with the predictions of various Monte 
Carlo generators and hadronization models. The dependence of the UE activity on 
the dimuon invariant mass is well described by PYTHIA and HERWIG++ tunes de- 
rived from the leading jet/ track approach, illustrating the universality of the UE ac- 
tivity. The UE activity is observed to be independent of the dimuon invariant mass in 
the region above 40 GeV/ c^, while a slow increase is observed with increasing trans- 
verse momentum of the dimuon system. The dependence of the UE activity on the 
transverse momentum of the dimuon system is accurately described by MadGraph, 
which simulates multiple hard emissions. 
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1 Introduction 

In hadron-hadron scattering, the "underlying event" (UE) is defined as any hadronic activity 
that cannot be attributed to the particles originating from the hard scattering, which is charac- 
terized by a large momentum transfer, or to the hadronization of initial- and final-state radi- 
ation. The UE activity is thus due to the hadronization of partonic constituents, not involved 
in the hard scattering, that have undergone multiple-parton interactions (MPIs) and to the 
hadronization of beam remnants that did not participate in other scatterings. These semihard 
interactions cannot be completely described by perturbative quantum chromodynamics (QCD) 
and require a phenomenological description involving parameters that must be tuned with the 
help of data [IJ. 

The experimental study of the UE probes various aspects of hadron production in high energy 
hadron-hadron collisions. In particular it is sensitive to the interplay of perturbative methods 
describing the hard process and phenomenological models of the soft interactions that attempt 
to simultaneously describe MPIs, initial- and final-state radiation, the colour flow between final 
state partons, and the hadronisation process. Understanding the UE in terms of particle and 
energy densities will lead to better modelling by Monte Carlo programs that are used in precise 
measurements of standard model processes and searches for new physics at high energies. The 
UE affects the estimation of the efficiency of isolation criteria applied to photons and charged 
leptons, and the energy scale in jet identification. It also affects the reconstruction efficiency 
for processes like H— > 77, where the primary vertex is partly determined from the charged 
particles originating from the UE. Hard MPIs are an important background for new physics 
searches, e.g. same-sign W production from MPIs [2J is a possible background to the same-sign 
double lepton SUSY searches |3|. 

The Compact Muon Solenoid (CMS) [4J, ATLAS, and ALICE experiments have carried out 
UE measurements at centre-of-mass energies (\/s) of 0.9 TeV and 7TeV using hadronic events 
(minimum-bias and single-jet triggered) containing a leading track-jet 13161 or a leading track IZl 
|8]. The analysis of the central charged particles and forward energy flow correlations in hard 
processes, e.g. pp — > W(Z)X — > £v{££)X [9J, provides supplementary insights into the nature 
of MPIs. In this paper, we use the Drell-Yan (DY) process llTOl with the muonic final state at a 
centre-of-mass energy of 7 TeV to perform a complementary UE measurement. The DY process 
with muonic final state is experimentally clean and theoretically well understood, allowing the 
particles from the UE to be reliably identified. The absence of QCD final-state radiation (FSR) 
permits a study of different kinematic regions with varying transverse momentum of 7*/Z 
due to harder or softer initial-state radiation (ISR). The comparison of the UE measurement in 
DY events with QCD events having a leading track-jet is useful for probing the UE activity in 
different processes. UE measurements using the DY process have been reported previously in 
proton-antiproton collisions at ^/s = 1.96 TeV ||TT| . 

The UE activity at a given centre-of-mass energy is expected to increase with the momentum 
transfer of the interaction. Events with a harder scale are expected to correspond, on average, 
to interactions with a smaller impact parameter and, in some models, to more MPIs |12t |13| . 
This increased activity is observed to reach a plateau for high energy scales corresponding to 
small impact parameter. In this paper we investigate some aspects of the UE modelling in 
detail by measuring the invariant mass dependence of the UE activity for DY events with small 
transverse momentum of the DY system. This measurement separates the scale dependence 
of the UE activity from the ISR effect. The universality of the model parameters, denoted as 
tunes, implemented in the various MC programs is tested by comparing their predictions with 
our measurements. The portability of the UE parameters across different event generators. 
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2 Observables 



combined in some cases with different parton distribution functions (PDFs), is investigated as 
well. The modelling of the ISR is studied by measuring the UE activity as a function of the 
transverse momentum of the DY system. Finally, the dependence of the UE activity on ISR and 
FSR is determined by comparing the measurements from DY events with previous results from 
hadronic events containing a leading jet where FSR also plays a role. 

The outline of the paper is as follows. Section 2 describes the various observables used in the 
present study. Section 3 summarizes the different MC models used and corresponding UE pa- 
rameters. Section 4 presents experimental details: a brief detector description, data samples, 
event and track selection criteria, correction procedure, and systematic uncertainties. Section 5 
presents the results on UE activity measured in DY events and the comparison with the mea- 
surements based on a leading track-jet. The main results are summarized in Section 6. 

2 Observables 

The UE activity is measured in terms of particle and energy densities. The particle density 
(l/[A?/A(A^)](Nch)) is computed as the average number of primary charged particles per 
unit pseudorapidity rj and per unit azimuthal separation Acp (in radians) between a track 
and the transverse momentum of the dimuon system. The pseudorapidity is defined as f/ = 
— ln(tan(0/2) ), where 9 is the polar angle measured with respect to the anticlockwise beam di- 
rection. The azimuthal angle (p is measured in the plane perpendicular to the beam axis. The en- 
ergy density (1/ [A?/A(A^)] (Spr)) is expressed in terms of the average of the scalar sum of the 
transverse momenta of primary charged particles per unit pseudorapidity per unit azimuthal 
separation. The ratio of the energy and particle densities, as well as the total charged-particle 
multiplicity N^h and the transverse momentum spectrum are also computed. The charged- 
particle multiplicity and transverse momentum distributions are normalized to unit area and 
to the average number of charged particles per event, respectively. Particles are considered 
as primary if they originate from the initial proton-proton interaction and are not the decay 
products of long-lived hadrons with a lifetime exceeding 10^^*^ s. Apart from the muons from 
the DY process, all charged particles in the central region of the detector with pseudorapidity 
< 2 and with transverse momentum pj > 0.5GeV/c are considered. 

The spatial distribution of the tracks is categorized by the azimuthal separation A^. Particle 
production in the away region (|A^| > 120°) is expected to be dominated by the hardest ISR 
emissions, which balance the dimuon system. The transverse region (60° < |A(^| < 120°) and 
towards region i\A(p\ < 60°) are more sensitive to soft emissions and, in particular, those due 
to MPIs. The relevant information about the hard and the soft processes is extracted from 
the tracking and the muon systems of the CMS detector and thus the derived observables are 
insensitive to the uncertainties of the calorimetric measurements. The DY events with dimuon 
mass Mpiii around the Z resonance are the least contaminated by background processes (heavy- 
quark, tt, W+jets, and DY — > tt production) ||T4l [151 and best suited for the measurement of 
the UE activity. 

The UE activity is studied as a function of the magnitude of the dimuon transverse momentum 
(Pt^ — I Pt 1 + Pt 2 1) ^ function of Mj,^,. The dependence of the UE activity on 

for high-mass dimuon pairs effectively probes the ISR spectrum. In order to minimize the 
background contamination, the p j^' dependence is studied only in the narrow mass window 
81 < Mppt < 101 GeV/c^. In contrast to the study of the UE activity in hadronic events using a 
leading track-jet IB HI, this energy scale is sufficiently large to saturate the MPI contributions. 
This observation is verified by studying the UE activity as a function of the dimuon mass in 
a wider mass range, where the total transverse momentum of the dimuon system is kept to a 
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minimum by requiring pj < 5 GeV/ c. 



3 Monte Carlo models 



The UE dynamics are studied through the comparison of the observables in data with various 
tunes of PYTHIA6 Ull and its successor PYTHIA8 (iZl UHl. MadGraph (version 5) flll |20| , 
which simulates up to six final-state fermions (including the muons), and POWHEG ||2T| . which 
includes next-to-leading-order corrections on the hardest emission, are also compared to our 
measurements. For these two generators, softer emissions are simulated by pj-ordered parton 
showers using PYTHIA6 tunes and matched with the hard process produced by the generators. 
Hadronization in PYTHIA6 and PYTHIAS is based on the Lund string fragmentation model [22]. 
The measurements are also compared to predictions of the HERWIG++ \23\ angular-ordered 
parton shower and cluster hadronization model [I24ll25|. 

The UE contributions from MPIs rely on modelling and tuning of the parameters in the MC 
generators. The MPI model of PYTHIA relies on two fundamental assumptions |IT2| : 

• The ratio of the 2— >2 partonic cross section, integrated above a transverse momen- 
tum cutoff scale, and the total of the hadronic cross section is a measure of the 
amount of MPIs. The cutoff scale por is introduced to regularize an otherwise di- 
verging partonic cross section. 



Here a/sq = 1.8 TeV and e is a parameter characterizing the energy dependence of 
the cutoff scale. 

• The number of MPIs in an event has a Poisson distribution with a mean that depends 
on the overlap of the matter distribution of the hadrons in impact-parameter space. 

The MPI model used here II26I includes showering of the MPI process, which is interleaved 
with the ISR. 

The tunes of the models vary mainly in the MPI regularization parameters, por and e, in the 
amount of colour reconnection, and in the PDF used. The Zl tune |127| of PYTHIA6 adopts the 
results of a global tuning performed by the ATLAS Collaboration II28II and uses the fragmen- 
tation and colour reconnection parameters of the ATLAS AMBTl tune |l29|. The parameters of 
the Zl tune related to the MPI regularization cutoff and its energy dependence are adjusted to 
describe previous CMS measurements of the UE activity in hadronic events f6\ and uses the 
CTEQ5L PDF The Z2 tune of PYTHIA6 is an update of the Zl tune using CTEQ6L1 [30], the 
default used in most CMS generators; the regularization cutoff value at the nominal energy of 
i/so = 1.8 TeV is optimized to 1.832 GeV/c. The value of the energy evolution parameter for the 
Z2 tune is 0.275, as for the Zl tune. The 4C [31 1 tune of PYTHIA8 follows a similar procedure 
as the ATLAS AMBTl tune, but includes ALICE multiplicity data as well. The values of the 
Pot(-\/so) ^ parameters for the 4C tune are 2.085 GeV/c and 0.19, respectively. The effective 
value of Pot at v^s = 7 TeV is about 2.7GeV/c for both the Z2 and 4C tunes. 



(Pt + Pot)^ ' 



(1) 



with 




(2) 
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4 Experimental methods 



The LHC-UE7-2 tune of HERWIG++ is based on ATLAS measurements of the UE activity in 
hadronic events [7j. The regularization cutoff parameter por for the LHC-UE7-2 tune is 3.36 GeV/c 
at a/s = 7TeV. The CTEQ6L1 PDF is used in conjunction with PYTHIA6 Z2, PYTHIA8 4C, Mad- 
Graph Z2, and HERWIG++ LHC-UE7-2, while CTIO t32J is used for POWHEG, and CTEQ5L for 
the PYTHIA6 Zl simulations. 

A comparison of these models with the measurements is presented in Section 5. 

4 Experimental methods 

The present analysis is performed with a sample of proton-proton collisions corresponding to 
an integrated luminosity of 2.2 fb^^, collected in March-August 2011 using the CMS detec- 
tor HH. 

Muons are measured in the pseudorapidity range < 2.4 with a detection system consist- 
ing of three subsystems: Drift Tubes, Cathode Strip Chambers, and Resistive Plate Chambers. 
Matching track segments from the muon detector to the tracks measured in the inner tracker re- 
sults in a transverse momentum resolution between 1% and 5% for pj values up to 1 TeV/c. The 
tracker subsystem consists of 1440 silicon-pixel and 15 148 silicon-strip detector modules, and 
it measures charged particle trajectories within the nominal pseudorapidity range \rj\ < 2.5. 
The tracker is designed to provide a transverse impact parameter resolution of about 100 fim 
and a transverse momentum resolution of about 0.7% for 1 GeV/c charged particles at normal 
incidence (rj = 0). 

The detector response is simulated in detail using the GEANT4 package Il33ll . The simulated 
signal and background events, including heavy-quark, tf, W+jets, and DY — > tt production, 
are processed and reconstructed in the same manner as collision data. 

4.1 Event and track selection 

The trigger requires the presence of at least two muon candidates. In periods of lower instanta- 
neous luminosity both muons were required to have pr > 7GeV/c, while in other periods the 
transverse momentum requirements were 13 GeV/c and 8 GeV/c for the leading and subleading 
muons, respectively. The trigger efficiency is above 95% for the offline selected DY events with 
the requirement of 81 < Mfifi < lOlGeV/c^. The offline selection requires exactly two muons 
reconstructed in the muon detector and the silicon tracker. Muon candidates are required to 
satisfy identification criteria based on the number of hits in the muon stations and tracker, 
transverse impact parameter with respect to the beam axis, and normalized X'^ of the global 
fit IITSl . The backgrounds from jets misidentified as muons and from semileptonic decays of 
heavy quarks are suppressed by applying an isolation condition on the muon candidates. The 
isolation variable I for muons is defined as 

I = {Z [pT(tracks) + £r(EM) + £r(HAD)] - 7r(AK)V} /p?^, (3) 

where the sum is defined in a cone of radius AR = {A(p)^ + {ArjY = 0.3 around the muon 
direction; Arj and A^ are the pseudorapidity and azimuthal separation between the muon and 
tracks or calorimetric towers. Here ^^(tracks) is the transverse momentum of tracks, excluding 
muons, with pj > \ GeV/c, £r(EM) is the transverse energy deposited in the electromagnetic 
calorimeter, £7 (HAD) is the transverse energy deposited in the hadronic calorimeter, and p 
is the average energy density l|34) in the calorimeter and tracker originating from additional 
inelastic pp interactions (pile-up) in the same bunch crossing as the DY interaction.The calcu- 
lation of p takes into account the number of reconstructed primary vertices in the event; the 
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average value of p is 5.6GeV/c. A muon is considered to be isolated if 7 < 0.15. Because of 
the energy density correction, the isolation efficiency is independent of the number of pile-up 
interactions. 

The selected muons are required to have opposite charges, transverse momenta larger than 
20GeV/c, and pseudorapidity \f]\ < 2.4. Both muons are required to be associated with the 
same vertex, which is designated as the signal vertex. The selected signal vertex is required to 
be within ±18 cm of the nominal interaction point as measured along the z direction. At least 
five tracks are required to be associated with the signal vertex, and the transverse displacement 
of the signal vertex from the beam axis is required to be less than 2 cm. These criteria select a 
pure sample of DY events with a total background contribution of less than 0.5% as estimated 
from simulated events. 

Tracks, excluding the selected muons, are considered for the UE measurement if they are well 
reconstructed in the silicon-pixel and the silicon-strip tracker, have pj > 0.5 GeV/c and \r]\ < 2, 
and originate from the signal vertex. To reduce the number of improperly reconstructed tracks, 
a high purity reconstruction algorithm [35J is used. The high purity algorithm requires stringent 
cuts on the number of hits, the normalized of the track fit, and the consistency of the track 
originating from a pixel vertex. To reduce the contamination of secondary tracks from decays 
of long-lived particles and photon conversions, the distances of closest approach between the 
track and the signal vertex in the transverse plane and in the longitudinal direction are required 
to be less than 3 times the respective uncertainties. Tracks with poorly measured momenta are 
removed by requiring (^(pr) / Pt < 5%, where (t{pt) is the uncertainty on the pj measurement. 
These selection criteria reject about 10% of primary tracks and 95% of misreconstructed and 
secondary tracks. The selected tracks have a contribution of about 2% from misreconstructed 
and secondary tracks. 

4.2 Corrections and systematic uncertainties 

The UE observables, discussed in Section 2, are corrected for detector effects and selection 
efficiencies. The measured observables are corrected to reflect the activity from all primary 
charged particles with transverse momentum pj > 0.5GeV/c and pseudorapidity \r]\ < 2. 
The particle and energy densities are corrected using a bin-by-bin technique. In the bin-by- 
bin technique, the correction factor is calculated by taking the bin-by-bin ratio of the particle 
level and detector level distributions for simulated events and then the measured quantity is 
multiplied by this correction factor. There is a small growth in the particle and energy densi- 
ties with increasing pj^ and in the towards and transverse regions. Because of this slow 
growth of densities the bin migration in pj^' and Mj,^, has a small effect on the measurements, 
therefore a bin-by-bin method is considered to be sufficiently precise. There is a fast rise in 
the energy and particle densities in the away region with the increase of pj^, but corrected re- 
sults using a bin-by-bin method are consistent with correction obtained from a Bayesian l|36| 
technique. The transverse momenta of the charged particles have very good resolution and are 
corrected using a bin-by-bin method. In this analysis the average of the calculated correction 
factors from PYTHIA6 Z2, PYTHIA6 D6T, and MadGraph Z2 is used to correct the experimen- 
tal distributions. The maximum deviation from the average correction factor is taken as the 
model-dependent systematic uncertainty, estimated to be 0.7-1.4% for the particle and energy 
densities. In the case of charged-particle multiplicity, there is substantial bin migration and the 
corrected results using the Bayesian [36 1 and bin-by-bin techniques differ by 10-15%. There- 
fore the charged-particle multiplicity is corrected using a Bayesian unfolding technique with 
a response matrix obtained using the PYTHIA6 Z2 tune. The systematic uncertainty related to 
the correction procedure is calculated by unfolding the data with response matrices obtained 
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4 Experimental methods 



using different tunes. 

In the analyzed data, there are on average 6-7 collisions in each bunch crossing. Tracks originat- 
ing from these pile-up interactions cause the UE activity to be overestimated, so the measure- 
ments are corrected for the presence of pile-up interactions. The correction factor is calculated 
as the ratio of the UE activity for simulated events with and without pile-up. The uncertainty 
in the modelling of the pile-up events is estimated by varying the mean of the expected number 
of pile-up events by ±1. This uncertainty in pile-up modelling affects the particle and energy 
densities by 0.3-1.0%. The effect due to pile-up events is small because only the tracks associ- 
ated with the same vertex as the muon pair are used. The results are also cross-checked with 
low pile-up 7TeV data collected during 2010 and the differences are found to be negligible. 

Table 1: Summary of the systematic uncertainties on the particle and energy densities (in percent). The 
first three rows show the systematic uncertainties for the particle density in the towards, transverse, 
and away regions. The last three rows report the systematic uncertainties for the energy density. The 
numbers outside the parentheses refer to the case where the densities are measured as a function of M^,^ 
and those in the parentheses correspond to the measurements as a function of p j''. 



Observable 


model 

M,,,, ip'j'') 


pile-up 

M„, ip'j'') 


isolation 

Mftf, (p't'') 


mis-ID 

M„, (p?.'') 


background 

M,.,, (pf ) 


total 

Mf,,. (p'7) 


l/[Ar]A{A^)]{Nci,) (towards) 
l/[Aj/A(Ai^)](Nch) (transverse) 
l/[A;?A(A^)](N,h) (away) 


0.8 (0.8) 
0.7 (0.9) 
0.7 (0.6) 


1.0 (0.9) 
0.9 (0.9) 
0.9 (0.3-0.9) 


0.9-1.5 (0.9) 
0.8-1.7 (0.8) 
0.8-1.6 (0.8) 


1.0 (1.0) 
0.9 (0.9) 
0.9 (0.9) 


0.7 (0.3) 
0.7 (0.5) 
0.5 (0.5) 


2.0-2.3 (1.8) 
1.8-2.3 (1.8) 
1.7-2.2 (1.5-1.7) 


1/[Aj/A(A(^)](Ept) (towards) 
1/[At]A{A(P)]{'Lpt} (trar\sverse) 
l/[Ar]A{Acj))]{T.pT) (away) 


1.2 (1.2) 
1.1 (1.4) 
1.0 (0.8) 


0.8 (0.7) 
0.7 (0.7) 
0.7 (0.3-0.7) 


1.1-2.0 (1.4) 

1.0- 2.5 (1.3) 

1.1- 2.2 (1.1) 


0.8 (0.8) 
0.8 (0.8) 
0.8 (0.7) 


0.8 (0.7) 
0.8 (0.9) 
0.7 (0.2) 


2.1-2.7 (2.2) 
2.0-3.0 (2.4) 
2.0-2.7 (1.6-1.7) 



We also consider possible systematic effects related to trigger requirements, different beam-axis 
positions in data and simulation, various track selection criteria, muon isolation, and misiden- 
tification of tracks. The combined systematic uncertainty related to trigger conditions, the vary- 
ing beam-axis position, and track selection is less than 0.5%. The systematic uncertainty due to 
isolation is calculated by removing the isolation condition in the simulated events used for the 
correction and is found to be 0.8-2.5% for the particle and energy densities. 

The yield of secondary tracks originating from the decay of long-lived particles is not correctly 
predicted by the simulation [37J. To estimate the effect of secondary tracks, a subset of simu- 
lated events is created by rejecting tracks that do not have a matching primary charged particle 
at the generator level. The uncertainty is evaluated by correcting the measurements with this 
subset of the simulated events, containing fewer secondary tracks, and is found to be 0.7-1.0% 
for the particle and energy densities. 

Though the total contribution of background processes is very small, it affects the measurement 
at higher (50-100 GeV/c) and small M^^ (40-60 GeV/c^) where the contamination from tt 
and DY— > rr background processes is 1% and 5%, respectively. The particle and energy den- 
sities differ between DY— )■ tt and DY— )■ fiji (the signal process) by 20%. The particle (energy) 
density for the tt background is two times (four times) that for the signal process. Combination 
of the differences in the densities for background processes and relative background contribu- 
tions gives a systematic uncertainty of 0.2-0.9%. 

Table [ijsummarizes the dominant systematic uncertainties on the particle and energy densities. 
The total systematic uncertainty on the particle and energy densities is in the range 1.5-3.0%, 
whereas the uncertainties on the track multiplicity and pr spectra reach 10% in the tail (not 
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reported in Table [Tjl. In all figures, inner error bars represent the statistical uncertainty only, 
while outer error bars account for the quadratic sum of statistical and systematic uncertainties. 

5 Results 

The UE activity in DY events, for charged particles with pj > 0.5GeV/c and \r]\ < 2.0, is 
presented as a function of M^^f, and pj^\ The multiplicity and the transverse momentum dis- 
tributions are also presented for two different sets of events, p^^ < 5 GeV/ c and 81 < M^,^ < 
101 GeV/c^. Finally, the UE activity in the transverse region is compared with that measured in 
hadronic events using a leading track-jet. 

5.1 Underlying event in the Drell-Yan process 

The energy-scale dependence of the MPI activity is studied by limiting the ISR. To accomplish 
this we require the muons to be back-to-back in the transverse plane with p^^' < 5GeV/c and 
measure the dependence of the UE activity on the dimuon mass, M^^,. The resulting particle 
and energy densities are shown in Fig. [ij Because the activity is almost identical in the towards 
and transverse regions, they are combined as | A^| < 120°. The contribution of ISR to the UE ac- 
tivity is small after requiring pj^' < 5 GeV/c, as shown by the prediction of HERWIG++ without 
MPIs. This figure also illustrates the dominant role of MPIs in our current models as they gen- 
erate more than 80% of the UE activity in these ISR-reduced events. The lack of dependence of 
the UE activity on Mj,^, within the range under study (40-140 GeV/c^) indicates that the activity 
due to MPIs is constant at energy scales down to 40 GeV. The quantitative description by model 
tunes based on the minimum-bias and UE observables in hadronic events is illustrated by the 
MC/Data ratios in Fig. [l] In general, PYTHIA6 Z2, PYTHIA8 4C, and HERWIG++ LHC-UE7-2 
describe the densities well, whereas the Z2 tune used together with the powheg generator un- 
derestimates both densities by 5-15%. Both PYTHIA and HERWIG++ model tunes derived from 
the UE measurement in hadronic events using the leading jet/ track approach describe the UE 
activity in the Drell-Yan events equally well and hence illustrate a certain universality of the 
underlying event across QCD and electroweak processes in hadronic collisions. 

Dependence of the UE activity on the transverse momentum of the dimuon system is shown 
in Fig. |2] in the towards, transverse, and away regions (top to bottom) for events having M^^, 
between 81GeV/c^ and lOlGeV/c^. At this high energy scale, the p^' dependence of the UE 
activity is sensitive to the ISR. The slope in the p^'^ dependence of the UE activity is identical 
for a model with and without MPIs and is therefore mainly due to ISR. The predictions of 
HERWIG++ without MPIs underestimate the measurements in the away region as well because 
the MPIs produce particles uniformly in all directions. The UE activity does not fall to zero 
when pj^ — > because of the presence of the hard scale set by M^^. 

The particle and energy densities in the away region rise sharply with p^^ and, because of 
momentum conservation mainly sensitive to the spectrum of the hardest emission, are equally 
well described by all tunes and generators considered. In the towards and transverse regions 
there is a slow growth in the particle and energy densities with increasing pj^'. The energy 
density increases more than the particle density, implying a continuous increase in the average 
transverse momentum of the charged particles with pj^. This effect is also reflected in the ratio 
of the energy density to the particle density. The activity in the towards region is qualitatively 
similar to that in the transverse region. Quantitatively, the activity is higher in the transverse 
region than the towards region, an effect caused by the spill-over contributions from the recoil 
activity in the away region, which balances the dimuon system. This observation is visible in 
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Figure 1: Top: The UE activity as a function of the dimuon invariant mass (M^^,) for events 
with < 5 GeV/ c for charged particles having Acp < 120° : (left) particle density; (centre) 
energy density; (right) ratio of the energy and particle densities. The predictions of PYTHIA6 
Z2, POWHEG Z2, PYTHIAS 4C, and HERWIG++ LHC-UE7-2 (with and without MPIs) are also 
displayed. In the top right plot, the structure around 60-80 GeV/c^ for HERWIG++ without 
MPIs reflects the influence of photon radiation by final-state muons, which is enhanced below 
the Z resonance. Bottom: Ratios of the predictions of various MC models and the measurement. 
The inner band shows the statistical uncertainity of data whereas the outer band represents the 
total uncertainty. 
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Fig. [2] at small pj^, where the radiation contribution is small and the activity in the transverse 
region is the same as that in the towards region. 

Figure |3] presents the ratios of the predictions of various MC models to the measurements for 
the observables shown in Fig. |2j Statistical fluctuations in the data induce correlated fluctua- 
tions for the various MC/data ratios. MadGraph in conjunction with PYTHIA6 tune Z2 de- 
scribes the dependence of the UE activity very well, both qualitatively and quantitatively. 
PYTHIA8 4C and HERWIG++ describe the p?^'' dependence of the particle density within 10-15%, 
but fail to describe the energy density. PYTHIA8 4C and HERWIG++ agree better with data as 
Pj^' approaches zero. The combination of the Z2 tune with powheg fails to describe the energy 
density in the towards and transverse regions, but gives a reasonable description of the particle 
density. This observation, combined with the information in Fig. [l| indicates that the discrep- 
ancies are not necessarily due to a flaw in the UE tune, but to an inadequate description of the 
multiple hard emissions and the different sets of PDFs used with POWHEG. At small p^^ the 
comparisons with PYTHIA6 Z2 and POWHEG Z2 are similar to those in Ref. ||38l, where PYTHIA6 
gives a good description of the pj^ spectrum while POWHEG underestimates the p^^. 

Figure |4] shows the distributions of charged particle multiplicity (top row) and transverse mo- 
mentum (bottom row). Figure |4] (left) shows a comparison of the normalized distributions 
in the away, transverse, and towards regions for events satisfying 81 < M^^, < lOlGeV/c^. 
As expected, the transverse and towards regions have fewer charged particles with a softer 
Pt spectrum than the away region. Figure |4] (centre) shows the comparison of the normal- 
ized distributions in the transverse region for two different subsets of the selected events, one 
with 81 < M,,,, < 101 GeV/c^ and one with pj^ < 5GeV/c. The charged particle multiplic- 
ity is decreased and the pr spectrum is softer when p^^ < 5GeV/c is required, because of the 
reduced contribution of ISR. Figure |4] (right) shows the comparison of the normalized distribu- 
tions with the predictions of various simulations in the transverse region for events satisfying 
81 < Mfipi < lOlGeV/c^. The charge multiplicity distribution is described well, within 10- 
15%, by MadGraph Z2 and pythia8 4C. The pr spectrum is described within 10-15% by 
MadGraph Z2, whereas pythia8 4C, powheg Z2, and herwig++ LHC-UE7-2 have softer 
Pt spectra. The various MC programs achieve a similar level of agreement with data in the 
towards region as in the transverse region. 

5.2 Comparison with the UE activity in hadronic events 

The UE activity was previously measured as a function of leading jet pj in hadronic events 
for charged particles with pseudorapidity \f]\ < 2 and with transverse momentum pj > 
0.5GeV/c [6|. Figure |5] shows the comparison of the UE activity measured in the hadronic 
and the DY events (around the Z peak) in the transverse region as a function of p'^^'^^^jet 
and pj^\ respectively. For the hadronic events two components are visible: a fast rise for 
pleading jet ^ iQGeV/c due to an increase in the MPI activity, followed by an almost constant 
particle density and a slow increase in the energy density with p^j^'^^^^ The increase in the 

particle and energy densities for p!^^*^™^ ^'^^ > 10 GeV/c is mainly due to the increase of ISR and 
FSR. Owing to the presence of a hard energy scale (81 < Mpipi < 101 GeV/c^), densities in the 
DY events do not show a sharply rising part, but only a slow growth with pj^ due to the ISR 
contribution. 

For pj^ and p!^^'^''^^ '"^^ > 10 GeV/c, DY events have a smaller particle density with a harder pj 
spectrum compared to the hadronic events, as can be seen in Fig.|5] This distinction is due to the 
different nature of radiation in the hadronic and DY events. Drell-Yan events have only initial- 
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Figure 2: The UE activity in the towards (upper row), transverse (centre row), and away (bot- 
tom row) regions as functions of p^'^ for events satisfying 81 < Mfifi < lOlGeV/c^: (left) particle 
density; (centre) energy density; (right) the ratio of the energy density and the particle densi- 
ties. Predictions of MadGraph Z2, powheg Z2, PYTHIA8 4C, and HERWIG++ LHC-UE7-2 
(with and without MPIs) are superimposed. 
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Figure 3: Ratios, as functions of pj^, of the predictions of various MC models to the measure- 
ments in the towards (upper row), transverse (centre row), and away (bottom row) regions for 
events satisfying 81 < M^^^ < 101 GeV/c^: (left) particle density; (centre) energy density; (right) 
the ratio of the energy density and particle densities. The inner band shows the statistical un- 
certainty on the data whereas the outer band represents the total uncertainty. 
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Figure 4: Distributions of the charged particle multiplicity (upper row) and transverse mo- 
mentum (bottom row) of the selected tracks. The left plots show the comparisons of the 
normalized distributions in the away, transverse, and towards regions for events satisfying 
81 < Mfiii < 101 GeV/c^. Comparisons of the normalized distributions in the transverse region 
are shown in the centre plots, requiring 81 < Mp,fi < 101 GeV/c^ or pj' < 5 GeV/c . The right 
plots show the comparisons of the normalized distributions in the transverse region with the 
predictions of various simulations for events satisfying 81 < Mj,^, < 101 GeV/c^. 



13 



state QCD radiation initiated by quarks, which fragment into a smaller number of hadrons 
carrying a larger fraction of the parent parton energy whereas the hadronic events have both 
initial- and final-state QCD radiation predominantly initiated by gluons with a softer fragmen- 
tation into hadrons. Similar behavior is observed for the track-jet measurement where the UE 
activity is higher by 10-20% for gluon-dominated processes, as estimated from simulation. 
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Figure 5: Comparison of the UE activity measured in hadronic and Drell-Yan events (around 
the Z resonance peak) as a function of p!^^'^™^^'^' and p^^', respectively: (left) particle density, 
(centre) energy density, and (right) ratio of energy and particle densities in the transverse re- 
gion. 



6 Summary 

We have used Drell-Yan events to measure the UE activity in proton-proton collisions at = 

7 TeV, which were recorded with the CMS detector at the LHC. The DY process provides a 
UE measurement where a clean separation of the hard interaction from the soft component is 
possible. After excluding the muons from the DY process, the towards (|A(/)| < 60°) and the 
transverse (60° < |A(^| < 120°) regions are both sensitive to initial-state radiation and multiple 
parton interactions. The DY process provides an effective way to study the dependence of the 
UE activity on the hard interaction scale, which is related to the invariant mass of the dimuon 
pair. The influence of the ISR is probed by the dependence on the transverse momentum of the 
muon pair. 

The UE activity is observed to be independent of the dimuon mass above 40 GeV/c^, after lim- 
iting the recoil activity, which confirms the MPI saturation at this scale. The UE activity in the 
DY events with no hard ISR is well described by PYTHIA6 and MadGraph with the Z2 tune 
and the CTEQ6L PDF. The Z2 tune does not agree with the data if used with PDFs other than 
CTEQ6L, as in the case of the POWHEG simulation. The PYTHIA8 4C and HERWIG++ LHC-UE7- 
2 tunes provide good descriptions of the energy-scale dependence of the UE activity. Thus 
the dependence of the UE activity on the energy scale is well described by tunes derived from 
hadronic events, illustrating the universality of MPIs in different processes. This universality 
is also indicated by the similarity between the UE activity in DY and hadronic events, although 
these events have different types of radiation. In addition, there is some ambiguity in the defi- 
nition of the hard scale for both types of events. 

The UE activity in the towards and transverse regions shows a slow growth with the transverse 
momentum of the muon pair and provides an important probe of the ISR. The leading-order 
matrix element generator MadGraph provides a good description of the UE dependence on 
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6 Summary 



dimuon transverse momentum. However, pythia, POWHEG, and HERWIG++, which do not 
simulate the multiple hard emissions with sufficient accuracy, underestimate the energy den- 
sity, but describe the particle density reasonably well. These measurements provide important 
input for further tuning or improvements of the Monte Carlo models and also for the imder- 
standing of the dynamics of QCD. 
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